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SUMMARY 


^ The  matrix  through- flow  analysis  is  extended  to  cover  turbulent  and  viscous 
flows  by  solving  the  streamwise  equation  of  motion  to  obtain  a suitable  entropy 
field  derived  from  the  stresses. 


The  results  of  computation  for  straight  circular  pipe  flows  agree  well 
with  published  data  for  laminar  flow  with  and  without  swirl  and  for  turbulent 
flow.  ^ 

V 

1 . Introduction 

Matrix  through-flow  analysis  is  widely  used  for  calculation  of  the  flow 
in  turbomachines.  It  solves  the  eqtiations  of  continuity,  motion,  energy  and 
state.  For  inviscid  fluids  these  equations  aire  sufficient  to  determine  the 
flow  pattern.  In  order  to  simulate  the  flow  of  a viscous  fluid,  a loss  model 
was  introduced  in  the  form  of  a local  polytropic  efficiency  by  Marsh^^.  The 
equations  remained  those  of  any  inviscid  flow. 

Marsh  and  Bosman^  derived  equations  which  are  consistent  with  the  loss 
model,  assuming  that  a dissipative  force  exists,  which  opposes  the  velocity 
vector.  The  loss  model  consists  of  either  empirical  data  or  equations 
relating  the  change  of  entropy  to  the  fluid,  the  flow  and  the  blade  geometry. 

The  present  paper  is  an  attempt  to  connect  the  changes  of  the  entropy  both 
to  the  viscous  terms  of  the  Navier-Stokes  equations  and  to  the  turbulence  terms 
of  the  mean  flow  Navier-Stokes  equations,  and  to  deduce  the  flow  of  a viscous 
fluid  without  including  empirical  equations  for  the  entropy.  It  is  based  on 
the  analysis  presented  by  Marsh  and  Bosman^ . 

2.  Flow  on  a Mean  Stream  Surface 


For  a three-dimensional  flow,  the  pattern  can  be  found  by  solving  the 
equations  of 

(a)  continuity 

(b)  motion  (3) 

( c ) energy 

(d)  state. 

niese  six  equations  are  sufficient  to  determine  the  fluid  properties  euid  the 
three  components  of  velocity. 

If  the  flow  takes  place  on  a prescribed  stream  surface,  then  there  is  a 
geometric  condition  relating  the  three  components  of  velocity,  which  must  be 
satisfied  in  order  that  the  velocity  vector  lies  on  the  stream  surface. 

This  geometric  condition  can  replace  one  of  the  equations  of  motion  in  the 
solution  for  the  flow  pattern. 


Replaces  A.R.C.37  01? 


The  replaced  equation  of  motion  can  be  used  to  evaluate  the  body  force 
which  must  be  applied  for  the  flow  to  follow  the  prescribed  stream  surface. 

If  a rectfmgulau?  co-ordinate  system  is  chosen  on  the  stream  surface  with 
axes  S along  Euid  N across  the  stream,  then  the  third  axis,  n,  is  normal  to 
the  surface  (Fig.l). 

The  body  force  F,  which  is  necessary  for  the  inviscid  flow  to  follow  the 
stream  surface,  is  normal  to  the  surface.  The  equations  of  motion  for  the  S 
and  N directions  therefore  do  not  contain  a component  of  F.  These  two 
equations  are  retained  for  calculating  the  flow  pattern.  The  equation  of 
motion  for  the  n direction,  which  is  replaced  by  the  geometrical  condition, 
contains  the  body  force  F and  it  can  be  used  to  determine  the  magnitude  of 
this  force. 

For  a reversible  adiabatic  flow  on  a prescribed  stream  surface,  the  flow 
pattern  is  determined  by  the  following  six  equations 

1.  The  continuity  equation. 

2a.  The  equation  of  motion  for  the  N direction  in  a co-ordinate  system 
rotating  with  the  blade  row 

N*W  X (Vx  V)  = N.(VI  - TVs  - D)  (l) 

or  since  N*D  = 0 

N«W  X (V  X V)  = N*(VI  - TVs)  (2) 

2b.  The  streamwise  equation  of  motion 

W*W  X (V  xV)  = W*(VI  - TVs  - D) 

and  using  the  energy  equation  for  steady  flow  we  get 
Ds 

T — = - W-D  (3) 

Dt 

2c.  A geometrical  condition  relating  the  three  components  of  velocity 
expressed  as 

W«F  = 0 (4) 


3.  Energy  equation  for  steady  flow 
DI 

— =0  (5) 

Dt 

assuming  that  it  is  valid  for  the  individual  streamlines  in  steady 
irreversible  adiabatic  flow. 

4.  The  streamwise  equation  of  motion  which  will  be  used  to  calculate 
the  entropy  field  from  the  dissipation  force  S.  This  is  given  by 

1 

I = - (V.  r) 

P 

Ibis  is  discussed  in  section  4. 


(6) 


This  approach  differs  from  the  one  by  Marsh  6md  Bosman  who  determine  the 
entropy  field  from  empirical  data  and  use  the  streamwise  equation  of  motion 
to  calculate  the  dissipation  force. 

3.  Analytical  Equations 

The  treatment  of  the  continuity  equation,  the  equation  of  motion  for  the 
N direction  and  the  geometric  condition  is  the  same  as  the  one  by  Marsh  and 
Bosman^  (where  the  bar  indicates  special  derivatives  , so  we  will  give  here 
only  the  final  equations. 


3.''  Continuity 

— fr  B p + — (r  B p = 0 

dr  \ J d%\ 

A stream  function  may  then  be  defined  as 


dif 

7-  = r B p 
dr 


dif 

— =-rB  p V 
ds  ' 


In  this  way  the  continuity  equation  is  satisfied. 
3.2  Geometric  conditions 
The  velocity  components  are  related  by 


W tam  X + W„  + Iff  tan  u = 0 
r u z 


where  by  definition 
tan  X = ^ 


tan  ^ = — 

”0 


□ 


3.3  The  N-direction  equation  of  motion 
This  finally  takes  the  form: 


<UIL  Ma.  it  jrLcmt 


lliese  are  defined  as 


0q 

dq 

n 

_ JL 

dq 

dq 

dq 

dq 

dr 

dr 

m 

XL 

dd 

ds 

dz 

dd 

(10a) 


(10b) 


and  they  represent  the  rate  of  cheuige  of  a quantity  q with  r or  z on  a stream 
surface  at  a given  value  of  z or  r respectively. 


I 


- If  - 


+ = V — (r  B p)  - V — (r  B p)  + (r  B p)®  — - 

ar®  as*  * ar  3z  dijf 

rBpfT  p 3s 


rBpT  r-  as  . , x"! 

- — (V_  - Wq  tan  ^l)  - — (V  - Wq  tan  X)  + 

W®  L ar  3a  -• 

r ? ^ \ n 

+ p B tan  — (r  Vq)  - tan  X — - (p  Vq). 

L ar  a®  J 


(11) 


3.4  The  S-direction  equation  of  motion 


'nils  is  given  as 


Ds 

T — = 
Dt 


W-D 


(4) 


where  D is  the  dissipation  force  and  it  is  opposite  to  the  velocity  vector, 
then 


D. 


■'e 


— = — = = K(W) 


(12) 


W 


w 


e 


w 


with  K<0  and  equation  (3)  becomes 


Ds 

T — = - W*K  (W). 

Dt 


(13) 


So  the  entropy  increases  along  a streamline. 

If  the  fluid  is  inviscid  the  dissipation  force  is  zero,  and  equation  (13) 
is  reduced  to  s = const  along  streamlines. 

4.  The  Dissipation  Force 

The  dissipation  force  as  shown  in  Appendix  1 is  given  as 

D = - (V.r).  (6) 

P 

This  can  be  calculated  for  every  point  in  the  field,  provided  that  the 
viscosity  is  known. 

Theoretically  equation  (12)  will  automatically  be  satisfied. 

In  practice  this  is  not  always  true,  because  of 

(a)  numerical  errors  involved  in  the  calculation  of  the  derivatives, 

(b)  approximations  involved  in  the  calculation  in  the  numerical  value 

given  for  the  viscosity. 


In  order  to  minimise  the  error  the  predominant  components  of  the 
dissipation  force  and  velocity  are  used  to  calculate  the  coefficient  K(w). 
This  is  expected  to  give  very  good  results  in  cases  where  there  is  a 
predominant  velocity  but  may  cause  problems  in  cases  without  it,  e.g. 
recirculating  flows. 


5-  The  Solution  of  the  Equations 


This  includes  the  following  steps: 


(a)  A distribution  of  the  stream  function  + is  assumed. 

(b)  From  equations  (8)  and  (9)  the  three  components  of  velocity  can 

be  calculated. 

(c)  Then  the  stress  components  can  be  calculated  from  equations  (12)  and 

(8)  respectively  of  Appendix  1. 


(d) 


The  streamwise 
entropy  field 
equation  (13) 


equation  of  motion  is  solved  in  order  to  get  the 
as  follows.  (See  Fig. 2).  Along  a streamline 
can  be  integrated, 


A 


Assuming  that  WK  is  constant  between  A and  B,  we  get 


Sg  = - W*K* (AB)/T. 


(14) 


can  be  calculated  as  a function  of  the  entropy  of  the 

neighbouring  points  and  (AB)  from  the  geometrical  definition  of 
the  stream  surface. 

(e)  The  right-hand  side  of  equation  (II)  cam  be  calculated  smd 

(f)  a bamd  matrix  may  be  formed  to  solve  equation  (II). 

(g)  The  resulting  + have  to  be  compared  with  the  initial  + , and  if 

the  difference  of  the  two  + is  greater  than  a certain  acceptable 
value  the  whole  procedure  has  to  be  repeated  from  the  beginning. 

6.  The  Flow  in  Axisymmetric  Pipes 

In  order  to  check  the  whole  procedure  a computer  program  was  written 
to  calculate  the  flow  in  a straight  circular  pipe. 

6.1  Equations  of  the  flow 

For  axisymmetric  flows 

- = 0.  (15) 

de 

"nie  special  derivatives  are  reduced  to  the  usual  partial  derivatives. 


Equation  (21)  can  be  integrated  along  a streamline.  The  angular 
momentiim  at  point  B (Fig. 2)  can  be  expressed  as 

r B K 

— dL 

(rVj  = (rv,)  . ■*  (22) 

B A 

Assuming  that  K is  constant  between  A and  B we  get 

VmJ 

(r  Vj  = (r  Vq)  e a (23) 

B A 

Equation  (23)  is  of  the  same  form  as  the  equation  derived  by  Talbot  (13)  for 
laminar  swirling  flows. 

6.2  The  computer  program 

This  follows  the  steps  described  in  section  5 with  an  extra  calculation 
for  the  solution  of  the  angular  momentum  equation  (e.g.  23)- 

6.2.1  The  finite  difference.-  Figs. 5 to  5 give  the  lattices 
used  for  the  finite  difference,  depending  on  the  position  of  the  point  relative 
to  the  boundaries. 

The  derivatives  are  replaced  with  the  following  expressions 


Coefficients  a.  . or  a„  are  calculated  using  Taylor's  series, 
ij  K, 

Uie  error  in  the  calculation  of  the  derivatives  is  of  the  order 

of— 1 

L dx*  J 

6.2.2  The  calculation  of  the  flow  conditions  upstream  of  a point. - 
In  order  to  solve  equations  (14)  and  (23)  the  conditions  at  the  point  A (Fig. 2) 
have  to  be  calculated.  Normally  point  A does  not  coincide  with  a grid  point 
so  some  approximations  had  to  be  made  for  the  determination  of 

(a)  entropy  at  A 

(b)  angular  momentum 
and  (c)  distance  AB. 

For  the  entropy  it  was  assumed  that 

S = S(i|r). 

A cubic-spline  polynomial  was  constructed  to  fit  the  entropy  at  i-1 
column  and  was  deduced  from  that  polynomial. 


- 8 - 


The  tangential  velocity  V-  was  assumed  to  vary  linearly  between  the 
points  which  include  A. 

The  distance  AB  was  approximated  by  a straight  line  between  A and  B. 

Hxe  radius  r^  was  found  again  by  cubic-spline  polynomial  fitting  between 
radius  r and  stream  function 

7*  Examples  of  the  use  of  the  Program 

7«1  Laminar  flow 

7-1 -1  Preservation  of  the  profiles.-  A fully  developed  two- 
dimensional  laminar  profile  was  set  as  the  initial  profile  at  the  entry  of 
the  pipe. 


7*'' ‘I  .I  Initial  conditions.-  Kie  velocity  is  given  by  the  relation 
\ (26.) 


V = 0. 
r 


(26b) 


Substituting  e.g.  (26a)  to  (l7a)  we  get 


— = 2 r p V 


(’-5) 


which  can  be  integrated  along  the  entrance  of  the  pipe  with  the  boundary 
condition 

+ = 0 at  r = 0 


to  give 


^ ° V 2 4a=  / 


Assuming  that  the  total  enthalpy  is  constant  at  the  entry, 


= 0 , 


we  get  equation  (l8) 


3 V r p T 3s 

p V + r p = p V — 

z 


— = . I ^ /Vn 

3r  T 3r  \ 2 / 


and  assuming  that 


s = 0 , 


V=2V  at  r=0 
z o 
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It  was  assmned  that  the  velocity  near  the  wall  has  a parabolic  distribution, 
as  a result  of  which  the  stress  component  is  given  by  the  relation 

'rz  = 2 pTz,  -VCr’  - .»). 

The  entropy  increase  along  the  wall  is  proportional  to  the  dissipation 
force  of  the  nearest  internal  point  or 

Sb  = S^-D^Wz)/T 


- 10  - 

7. 1.1 .4.  The  results.-  Fig. 6 gives  the  calculated  profiles  and  the 
analytical  ones. 

For  a distance  of  20  diameters  the  velocity  profile  changes  by  about 
2.551^  of  the  analytical  value. 

lliis  deviation  can  be  attributed  to  the  fact  that  the  order  of  magnitude 

r n 

of  the  error  of  the  finite  difference  approximation  is  0 — and 

L J 


this  is  not  zero  as  can  be  seen  from  equation  (28). 


7.1.2  The  development  of  laminar  flow.-  An  initial  uniform  velocity 
profile  was  set  at  the  entry  of  the  pipe  and  the  development  of  the  flow  profile 
was  predicted. 

7. 1.2.1  Initial  conditions.-  The  distribution  ci  velocity  across 
the  entrance  is 


V = V 
z 0 


(32a) 


V = 0 

r 


(32b) 


Substituting  equation  (32a)  into  (I7a)  we  get 

(V 

— = r p V 

ar 

The  above  equation  can  be  integrated  along  the  entrance  of  the  pipe  with  the 
boundary  condition  + = 0 at  r = 0 


to  give 


♦ = p^o7 


Across  the  entrance  of  the  pipe  we  again  take 

H = const, 
o 

and  for  the  entropy  | 

S(r,0)  = 0. 

7. 1.2. 2 Boundary  conditions.-  Table  2 gives  the  boundary  conditions 

applied. 

It  is  assumed  again  that  the  velocity  near  the  wall  has  a parabolic 
distribution. 


rf  — — 
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Table  2 

Boundary  conditions  for  the 
development  of  laminar  flow 


Quantity 

Value  on 
the  axis 

Value  at 
the  wall 

Stream  function 

♦ 

+ = const. 

Axial  velocity 

V 

z 

Hi.  0 

Sr 

0 

Radial  velocity 

V 

r 

0 

0 

^rz 

0 

2 Vzy  Va/(r*  - a*) 

’^zz 

sn 

0 

Increase  of 
entropy 

0 

^ Vz  ’’rzw 

At  a 

7.1 .2.4  The  results.-  Fig. 7 gives  the  calculated  axial  velocity 
profiles  at  different  distances  from  the  entry,  with  the  corresponding 
experimental  profiles  by  Nikuradze  in  Schlichting^^, 

There  is  a discrepancy  of  about  (>%  for  the  values  of  the  velocity  at  a 
distance  of  5 diameters  which  is  reduced  gradually  to  almost  zero  at  distances 
higher  than  25  diameters.  We  think  that  this  is  due  to  the  inaccuracy  of 
calculating  the  shear  stresses  near  the  wall  at  the  first  few  stations  downstream 
of  the  entrance  of  the  pipe. 

Nevertheless,  the  values  of  the  velocity  at  the  centreline  are  identical 
with  those  calculated  by  Boussinesq  as  they  are  given  by  Schlichting^^  on  the 
graph  with  Nikuradze 's  date. 

7.1.3  Effect  of  swirl  on  fully  developed  laminar  flow.-  A 
tangential  velocity  is  superimposed  on  a fully  developed  laminar  flow  at  the 
entry  of  a pipe,  and  the  development  of  the  flow  profile  is  given. 


12  - 


relation 


7.1 O'l  Initial  conditions.-  The  axial  velocity  is  given  by  the 


V = 2 V 

2 O 


(35a) 


(35b) 


(35c) 


The  stream  function  distribution  is  given  by  the  relation 


♦ = - - — 


The  total  enthalpy  is 


= const. 


and  the  distribution  of  entropy  is 


1 r»  / r*  N 

j(r,0)  = - 4 V » — ( 1 ) - 

T L a**  \ 2a*  / 


7. 1.3. 2 Boundary  and  exit  conditions.-  These  are  the  same  as  for 
the  fully  developed  laminar  flow,  except  that  at  the  axis  we  now  have 


or  the  axial  velocity  at  the  centreline  is  equal  to  the  velocity  at  the  nearest 
point  with  the  same  z. 

7. 1.3. 3 The  results.-  Fig. 8 gives  the  axial  velocity  profile  for 

an  initial  swirl  00-  4 at  a distance  of  5 diameters  from  the  entry  of  the  pipe. 

At  the  axis  the  velocity  is  reduced  and  the  velocity  towards  the  pipe 
wall  is  increased. 

13 

Similar  behaviour  is  shown  theoretically  by  Talbot  , although  he  predicts 
smaller  values  for  the  velocity  at  the  centreline. 

Fig. 9 gives  the  flow  profiles  for  a distance  of  25  diameters  and  for  (0=  4. 

The  axial  velocity  profile  becomes  flatter  at  the  centreline  as  the 
angular  velocity  increases. 

For  values  of  oaa/Vo  > 1 for  Re  between  20  and  2000  the  computer  program 
breaks  down.  This  is  due  to  the  fact  that  during  the  computation  reverse 
flow  appears  at  the  axis  which  moves  towards  the  entry  from  iteration  to  iteration. 
The  program  breaks  when  the  reverse  region  reaches  the  entry,  because  a 
discontinuity  of  velocity  appears  at  the  entry  with  which  the  initial  conditions 
cannot  cope. 


f 
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Fig. 10  shows  the  decay  of  the  tangential  velocity  profile.  The  profiles 
are  similar  for  different  o),  or  the  ratio 

V“. 

is  independent  of  the  angular  velocity  <0. 

The  profiles  are  similar  to  those  given  by  Lavan  et  al^. 

Fig. 11  gives  the  distribution  of  the  radial  velocity  at  a distance  of 
8 diameters  from  the  entry,  for  3 angular  velocities.  The  profiles  are  similar 
to  those  given  by  Talbot ^5. 

7.2  Turbulent  flow 

The  equations  of  motion  as  described  in  vector  form  by  equation  (1)  or 
in  analytical  form  in  Appendix  1,  remain  the  same  for  the  mean  velocities. 

The  only  difference  is  that  the  dissipation  force  now  contains  derivatives  of 
the  Reynolds  stresses  due  to  turbulence,  and  an  extra  problem  appears,  namely, 
how  to  calculate  these  extra  stresses. 

One  basic  assumption  made  by  those  who  have  used  numerical  models  for 
turbulence  is  that  the  turbulence  is  isotropic.  Hence  it  follows  that  there 
is  a scalar  turbulent  viscosity  which  satisfies  the  relation 

r 3 U . 3 U i -I 

i ^ (36) 

3xj  dx^  J 


- o u.  u.  = 


and  also  links  the  stresses  due  to  turbulence  with  the  mean  velocity  gradients 
of  the  flow  at  each  point. 


For  the  pipe  flow  an  initial  mean  velocity  distribution  is  assumed 

Vb 


= 


■’m 


(^) 


(37) 


The  distribution  of  the  stream  function  across  the  entrance  of  the  pipe 
will  be  IB/  iB/g 

r—  a 

V, 


♦ = p 


9V 

a 


p (a  - r)  a 

^ - -(.-r) 


9^ 

190  J 
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the  boundary  condition  for  ♦ being  that 

♦ = 0 at  r = 0 and  z = 0 . 

The  initial  distribution  of  entropy  is  given  by  equation  (31) 


(38) 


(59) 


The  boundary  conditions  and  the  exit  conditions  are  the  same  as  for  the 
fully  developed  laminar  flow,  with  the  exception  of  the  stress  t at  the  wall. 


r 
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|i  = c a u 


7.2.1  Solution  with  given  As  a first  step  towards  a solution 

of  the  turbulent  flow  in  a pipe,  the  distribution  of  the  turbulent  viscosity  was 
taken  as 

(40) 

(4la) 
(4lb) 


where  c takes  the  following  values 

c = 0.07  (a  - r)/a  for 

auid  c = 0.07  for 

u*  = 0.035 


0 < a - r < 0.25  a 
0.25  a < a - r < a 


Q 

These  were  taken  from  Hinze  based  on  data  by  Laufer  for  Re  = 50  000. 
The  wall  shear  stress  is 

. *a 


rz 


= p u 


(42) 


7. 2. 1.1  Results.-  Fig. 12  gives  the  profile  of  the  axial  velocity 
across  the  entrance  of  the  pipe  and  also  after  a distance  of  6 diameters. 

The  maximum  deviation  is  of  the  order  of  0.4^. 

7.2.2  The  K -e  model  for  turbulence.-  A more  elaborate  model  for 
turbulence  is  the  K -e  model  proposed  first  by  Harlow  and  Nakayama°. 

Hanjalic^  developed  two  transport  equations,  one  for  the  kinetic  energy  of 
turbulence  K,  and  one  for  the  rate  of  dissipation  of  turbulence. 

In  cylindrical  co-ordinates  for  axisymmetric  flow  and  for  high  Reynolds 
numbers  the  equations  are; 


De 

1 

- 1 d 

d 

de  \- 

e 

Dt 

P 

_ r dr 

dr  ) 

1 

dz 

dz  /_ 

+ 

Cl  — 

P 

® 

K 

D/C 

1 
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d 

dK  N 
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e 

Dt 

P 
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1 
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P 
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\\  i 

•1 

^ ( 

12^ 

12r 

Cr  s 

2 

V dr  J 

J 

* \ 

dr 

T 

dz  . 

J \ 

\ * ( 

♦ 

\ dr  J 

dz  J 

K 


(45b) 


(7)* 


(44) 


and  the  turbulence  viscosity  will  be 

= Ou  P 


(45) 
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nie  various  constants  appearing  in  these  equations,  as  given  by  Launder  and 
Spalding^,  are 

= 0.09,  = 1.44,  Cjj  = 1.92 

<r^  = 1 , <r  = 1.3  , *^1)  = 0.08. 

Equations  (43a)  and  (43b)  have  to  be  solved  for  each  iteration  of  the  stream 
fxinction. 


The  initial  conditions  are 


K 


e = K ft 

taken  from  Launder  and  Spalding^ . 

For  points  near  the  wall  it  is  assumed  that  the  velocity  profile  follows 
the  logarithmic  law. 


'zw 


u„ 


£r\ 


1^9  u^(a 


‘T  0.41 
Then  the  kinetic  energy  of  turbulence  will  be 

Kw  = u.^0.3 

and  the  rate  of  dissipation  will  be 


e,  = /C‘-»  /[0.41  (a-  r^  ] . 

7.2.2. 1 Results.-  The  initial  mean  velocity  profile  does  not 
change  along  the  pipe  by  more  than  -O.Jd(>  for  a distance  of  6 diameters  downstream 
of  the  entrance. 


8.  Conclusions 


A new  approach  to  the  use  of  the  dissipation  force  in  through-flow 
analysis  is  discussed  which  permits  more  accurate  calculation  of  a flow  field 
when  no  data  for  the  entropy  exists. 

The  dissipation  force  is  calculated  taking  into  account  the  viscous  £ind 
turbulence  terms  of  the  Navier-Stokes  equations. 

A computer  program  was  written  and  different  cases  of  pipe  flow  were 
checked . 

The  fully  developed  laminar  flow  profile  remains  unchanged  to  within 
2.5^  for  a distance  of  20  diameters. 


] i 


1 


Hie  development  of  the  laminar  flow  is  predicted  satisfactorily. 


Hie  effect  of  swirl  on  fully  developed  laminar  flows  is  predicted,  but 
because  of  lack  of  suitable  experimental  data  it  has  not  been  possible  to 
compare  the  results  with  experiments. 

In  the  case  of  turbulent  flow,  two  models  for  the  Reynolds  stresses  were 
used,  both  with  very  satisfactory  results. 
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Nomenclature 


B 

D 

F 

h 

H 

o 

I 

K 

K 

n 

N 

r 

s 

S 

T 

V 
W 
z 

e 

(d 

T 

L 

£ 

V 


integration  factor  for  the  continuity  equation 
dissipation  force  per  unit  mass 
body  force  per  unit  mass 
enthalpy 

stagnation  enthalpy 

rothalpy,  I = - o)  rVg 

ratio  of  dissipation  force  to  velocity 

kinetic  energy  of  turbulence 

normal  to  stream  surface 

normal  to  streamline 

radius 

entropy 

streamwise  direction 
temperature 
velocity  vector 
relative  velocity  vector 
axial  direction 
circumferential  direction 
angles  defined  by  equations  (10) 
angular  velocity 
stress  tensor 

distance  along  a streamline 

length  scale  of  turbulence 
meridional  velocity 


I 

I 

i 
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♦ 

a 

V 

0 
a 

V 

i> 

Re 

a 

X. 

1 

Subscripts 

w 

CL 

i,  j,k 


stream  function 

coefficients  for  the  calculation  of  the  derivatives 
mean  velocity  of  the  laminar  flow 
radius  of  the  pipe 
kinematic  viscosity 
a representative  variable 

dVo 

Reynolds  number,  Re  = 

p 

maximum  velocity  at  the  axis  for  turbulent  flow 
tvirbulent  viscosity 

rate  of  dissipation  of  tiurbulence  energy 
Cartesian  space  co-ordinate 


point  near  the  wall 
point  on  the  axis 

subscripts  denoting  Cartesian  co-ordinate  directions. 
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Appendix  1 

The  Caloulatlon  of  the  Dissipation  Force 
1 . The  Havler-Stokes  Bauatlone 

The  momentum  equations  in  cylindrical  co-ordinates  are  given  iDy  Moulden 
and  Wu^^  for  steady  flows 
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00 
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P 

0Z 

- - - (V-r), 


— -(V-r), 


= - - — + - (Vr), 


(iv.1  a) 

(A1b) 

(a1c) 


Using  the  relations 


I = H^-wrV.  = h + ^V>-wrV, 


(A2) 


and 


we  get 


7h  = T 7 s + — 7p 
P 


The  equations  of  motion  as  used  by  Bosman  and  Marsh  are 


— - T — - D 
dz  dz  ^ 


Comparing  equations  (A4)  and  (A5)  we  get 


fe 


2,  The  Components  of  the  Stress  Tensor 

k j kk ^ ^ 0 j-i 


According  to  Moulden  and  Wu  the  expression  (7»t)  is  given  hy  the 
following  relations: 


d ^ d d 

(V.t)^  = ““  + - — T « + — T + 

r . zr  30^0  - rr 

dz  r do  dr 


^rr  " ’’ee 


1 d 


(V*T)n  = T-.A  + - ” + 2 


r d0 


a d T , 

(Vt).  = — t „ + ■ ’’.A  + — + 


dz  r de  dr 

The  strain  tensor  is  given  hy  the  relations 

d V_ 


“ ’'rz  

dr  r 


1 3 Wq 


r d0 


Lls. 

dz 

3 1^9 

1 d V 

z 

dz 

+ • ■“ 

r d0 

Lh 

dVr 

dr 

dz 

1 LIr  ^ \ 


r d0 


For  fluids  which  follow  Stokes*  law  the  stresses  are  given  ty  the 
following  relations 


r = - ^ ^ (V.V)  + ji  e 

rr  rr 


’■flfi  = - I M e. 


’■*7  * " § ■*“ 
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II 

^zr 

’■*0  = 

’’da  = 

»^®0z 

II 

’"er  = 

^^®0r 

where  the  divergence  of  the 

velocity 

vector 

7*^  = 

d V 

z 

1 

+ 

dz 

r 30 

• For  incompressible  fluids  the  divergence 
equal  to  zero  from  the  continuity  equation. 

So  if  p = const. 

V*V  = 0 
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Lattice  points 
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Fig.  3*  Main  Lattice. 


LL 


I.L-U. 


I.U 


1- 


Fig.  4.  Lattice  for  points  near  the  weG.1. 


Fig.  5*  Lattice  at  the  entry. 
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g.  6.  Laminar  flow  profile 


Fig.  7fe»  The  development  of  the  pipe  flow 


Fig.  7c.  The  development  of  the  pipe  flow 
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